Self-consistent slave rotor mean field theory for strongly correlated systems 
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Building on work by Florens and Georges, we formulate and study a self-consistent slave rotor 
mean field theory for strongly correlated systems. This approach views the electron, in the strong 
correlation regime, as a composite of a neutral spinon and a charged rotor field. We solve the coupled 
spinon-rotor model self-consistently using a cluster mean field theory for the rotors and various 
ansatzes for the spinon ground state. We illustrate this approach with a number of examples relevant 
to ongoing experiments in strongly correlated electronic systems such as: (i) the phase diagram of 
the isotropic triangular lattice organic Mott insulators, (ii) quasiparticle excitations and tunneling 
asymmetry in the weakly doped cuprate superconductors, and (iii) the cyclotron mass of carriers in 
commensurate spin-density wave and U(l) staggered flux (or d-density wave) normal states of the 
underdoped cuprates. We compare the estimated cyclotron mass with results from recent quantum 
oscillation experiments on ortho-II YBCO by N. Doiron-Leyraud et al (Nature 447, 565 [2007]) 
which appear to find hole pockets in the magnetic field induced normal state. We comment on the 
relation of this normal ground state to Fermi arcs seen in photoemission experiments above Tc. 
This slave rotor mean field theory can be generalized to study inhomogeneous states and strongly 
interacting models relevant to ultracold atoms in optical lattices. 
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I. INTRODUCTION 

It is well known that strong electronic correlations in 
systems such as the transition metal oxides and ultra- 
cold atoms in optical lattices can lead to a wealth of 
novel phases and interesting phenomena not observed in 
conventional metals and superconductors. A continuing 
theoretical challenge is to develop tools to solve model 
Hamiltonians, such as the Hubbard model or the t-J 
model, to gain a qualitative as well as quantitative under- 
standing of the strong correlation regime of these models 
where perturbative methods are no longer expected to 
work. A variety of different approaches have been de- 
veloped to address this issue, and they can be classified 
roughly into two categories. The first category includes 
numerical methods such as exact diagonalization which 
are good at capturing the short distance physics^. A 
much more sophisticated approach is dynamical mean 
field theory (DMFT) which maps the Hubbard Hamil- 
tonian onto a quantum impurity model self-consistently 
coupled to a dynamical bath^i^. This approach and its 
extensions, which include some degree of spatial fluc- 
tuations, have led to significant progress in our under- 
standing of the Mott (metal-insulator) transition as well 
as the short-range spin correlations and superconduc- 
tivity in doped Mott insulators^. In the second cate- 
gory are semi-analytical approaches including renormal- 
ized mean field theorj^j^iS' and slave particle mean field 
theories. Slave-particle mean field theories, such as slave 
boson^iiiSiii and slave rotor~ theories, are motivated by 
the observation that strongly interacting systems often 
display (i) quite disparate charge and spin response with 
(ii) little evidence of sharply defined electron-like quasi- 
particles. This observation has led to the phenomeno- 
logical spin-charge separation postulate that electrons in 
many strongly correlated materials may be better viewed 
as composites of 'chargons' and 'spinous', which respec- 
tively describe the charge and spin degrees of freedom. 



A physically transparent way to encapsulate this idea is 
provided by the slave particle approaches which enlarge 
the Hilbert space of the physical electronic degrees of 
freedom in terms of separate chargon and spinon Hilbert 
space. The unphysical states are then eliminated by en- 
forcing constraints on the enlarged Hilbert space. The 
strongly correlated model then maps onto a model of 
interacting slave particles coupled to gauge field o^^'^'^ . 
This rewriting of the strongly correlated electronic model 
is useful in situations where the gauge fluctuations are 
weak. This assumption, that the gauge fluctuations are 
weak and can be ignored, is a reasonable approximation 
for conducting (metallic or superconducting) states, bro- 
ken symmetry insulating phases such as the Neel anti- 
ferromagnet, or if the gauge fluctuations are gappecU^. 
In other cases, we have to appeal to comparisons with 
experiments to see if such slave particle mean field ap- 
proaches provide a more useful starting point to think 
about the physics compared to perturbative analyses. 

In this paper we present a further development in 
the slave particle mean field approach. In conventional 
U(I) slave boson theory for the ground state of the t- 
J modeliSiiii, one is in the regime where the on-site 
Coulomb repulsion U ^ t, and the bosons are assumed 
to be condensed with the boson fluctuations being ig- 
nored. However, the boson model is a strongly inter- 
acting model which has to be solved in order to obtain 
its ground state and excitations. Similarly, for systems 
close to a Mott transition, number (charge) fluctuations 
are signiflcant so it is more appropriate to describe the 
charge degree freedom using slave rotorsi^. The rotor 
model is also nontrivial: its Hamiltonian resembles the 
Bose-Hubbard model and again it has to be solved. For 
instance, at commensurate densities, the rotor model can 
undergo a Mott transition leading to an electronic Mott 
transition in the coupled spinon-rotor problem. We study 
a number of simple situations and present results which 
suggest that solving the rotor model and coupling it self- 
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consistently to the spinons provides a reasonable mean 
field framework to study strongly correlated electronic 
systems, and also recovers U(l) slave boson mean field 
results for [//< ^ oo. 

The rotor sector have been studied by several authors 
using single site mean field theory as originally applied 
to the Bose-Hubbard modeli^ii^. We shall show that 
while the single site approximation captures the essence 
of Mott transitioiii^, it is inadequate in some respects and 
sometimes yields incorrect results. Going beyond this re- 
quires the evaluation of short range spatial correlations. 
We propose to solve the rotor Hamiltonian by diagonaliz- 
ing it on a small cluster coupled to the rest of the lattice 
self-consistently through an order parameter bath. We 
demonstrate that such a cluster mean field theory can get 
rid of the aforementioned problems of single site approx- 
imation while keeping the whole procedure conceptually 
simple and computationally affordable. This constitutes 
the main new aspect of the present work. This procedure 
is rather general, and we expect it would also be useful 
for other strongly interacting bosonic models, e.g., those 
which arise in the context of cold atoms in optical lat- 
tices. It can also be adapted to address inhomogeneous 
electronic states. In this paper, we apply this technique 
to address issues arising from some recent experiments in 
strongly correlated electronic systems. 

The outline of this paper and a summary of the main 
results is as follows. We begin in sections II with a re- 
view of the basic idea of the slave rotor approach and 
the tU J model Hamiltonian. In section III we introduce 
the self-consistent mean field theory of the rotor sector 
and spinon sector. Sections IV and V discuss various 
mean field ansatzes for the spinon Hamiltonian, and the 
cluster mean field theory for the rotors. The remain- 
der of the paper deals with applications of the coupled 
spinon-rotor mean field theory. In Section VI, we study 
a metal to Mott insulator transition as well as the possi- 
bility of a spin liquid phase in a Hubbard-type model at 
half-filling on triangular lattice. This model describes or- 
ganic compounds, such as k-(ET)2Cu2(CN)3, which fea- 
ture a pressure-tuned Mott transitio n^^i^^ . We show that 
in contrast to a single-site mean field theory, the cluster 
approach leads to the possibility of a spin liquid state 
between the metallic state and the antiferromagnetically 
ordered insulator. Section VII illustrates the apphcation 
of the mean field theory to the superconducting state 
of doped cuprates. Specifically, we find that properties 
of quasiparticle excitations such as their spectral weight 
and Fermi velocity, which are measured in photoemis- 
sion experiments, are in qualitative agreement with vari- 
ational wavefunction and renormalized mean field theory 
results. We also study the particle-hole asymmetry ob- 
served in recent tunneling experiments^^ on the weakly 
doped cuprates, and make an estimate of how this asym- 
metry varies with the cutoff energy scale in tunneling 
asymmetry imaging. In section VIII, we study a strongly 
correlated commensurate spin-density wave (SDW) state 
and a U(l) staggered flux or d-density wave (DDW) state 



as possible underlying normal states of the underdoped 
cuprates. We use the slave rotor mean field theory to 
estimate the effective cyclotron mass of the carriers in 
these states. These results are compared with recent high 
field quantum oscillation experiments on ortho-II YBCO 
which shows evidence of hole pockets. We suggest that 
doping dependent studies of the cyclotron mass should 
be able to distinguish between these two normal states. 



II. MODEL AND SLAVE ROTOR 
REPRESENTATION 

We consider the tUJ model for electrons on a two di- 
mensional lattice 

H = - ^v4aCja + ^ H '^iT"a + J^^^■^j (1) 

Here = t, t' represents hopping amplitudes to near- 
est and next nearest neighbor sites respectively. U > 
is the local Coulomb repulsion. Si = c\^da()C^ij is the 
electron spin operator, where a = {<Jx,o'y,o'z) are the 
Pauli matrices. J > is an antiferromagnetic spin ex- 
change between spins on neighboring sites. Ordinarily, 
the superexchange J derives from a combination of the 
kinetic terms and the Hubbard repulsion in the regime 
U ^ t,t', with J ^ At'^/U. Here, however, we retain 
an explicit J term in the model since we will analyze 
the spin physics of this model below using mean field 
theory. The J term is required in order to reproduce 
the strong coupling antiferromagnetic spin correlations 
within mean field theory (which decouples the spin and 
charge degrees of freedom). Even for weak to moderate 
couphngs, U < t, the Hubbard repulsion leads to an- 
tiferromagnetic spin fluctuations: in this regime, the J 
term can be viewed as a way to reproduce the weak cou- 
pling equal time spin structure factor within mean field 
theory. The tUJ model have been studied previously by 
several authorsi^ii^SiSi. We will solve this model using the 
slave rotor method, for a range of densities around half- 
filling and for a range of repulsive interactions U. For 
studying interaction driven transitions such as the Mott 
transition, we will keep J fixed for simplicity. In order to 
study the doping dependent correlations at fixed U, we 
set the antiferromagnetic exchange J — At^/U. 

We begin with a review of the slave-rotor description of 
strongly correlated electron systems introduced by Flo- 
rens and Georgesi^. The electron Hilbert space at a sin- 
gle lattice site has four states, |0), | t), | i), | Ti)- That is, 
it could be empty, or be occupied by a spin up electron, or 
a spin down electron, or be doubly occupied. In the slave 
rotor representation, the electron charge degree of free- 
dom is described by a charged rotor and the spin degree 
of freedom is described by a spin-1/2 spinon (fermion). 
Each of the four physical states then corresponds to a 
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direct product of the rotor state and the spinon state, 
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mean field theory as well as a sigma model for the rotor 
fluctuations) and the Coulomb blockade phenomenon in 
quantum dots^^. Lee and Lee have used the slave rotor 
representation to study possible spin liquid phases near 
the Mott transition in Hubbard model on the triangular 
lattice^^. 



Here on the r.h.s., the first ket \-nP) is the eigenstate of 
rotor charge with eigenvalue nP = 0,±1, and the second 
ket is eigenstate of spinon occupation number, nf^cr = 0, 1 
for a =ti i- Notice we have chosen a background charge 
+1 for the state with no electrons and each added elec- 
tron contributes charge —1. The enlarged rotor-spinon 
Hilbert space contains unphysical states such as |1)| t)- 
These unphysical states are avoided by imposing the op- 
erator constraint 



III. SLAVE ROTOR MEAN FIELD THEORY 

In the spirit of slave particle mean field theory, we ap- 
proximate the electron wavefunction by the direct prod- 
uct of appropriate spinon and rotor wavefunctions, with 
the constraint equation satisfied on average in the result- 



ing mean field state. 



1. Together 



with 



= (n^), this leads to two mean field constraints, 



(2) 



In the slave rotor representation, the electron number is 
equal to the spinon number, i.e., 

<^ = (3) 
The electron creation (annihilation) operator 

4. = 4.^"^'-' (4) 

c.,. = /..e+^'^N (5) 



where fa is the spinon annihilation operator, and the 
rotor creation (annihilation) operator , e"'"**' (e~**'), is 
defined by 



e±'^Mn?) = |n?±l) 



(6) 



We can then rewrite the tU J Hamiltonian in terms of the 
spinon and rotor field operator as 



^-Y^nt{nt-1) + JY^^{-^. (7) 



Here we have expressed the Hubbard repulsion between 
electron charges in terms of the rotors since only the ro- 
tors carry charge. Similarly, the antiferromagnetic ex- 
change interaction is expressed in terms of the spinous, 
with S{ = fj^ffa/sfifj since only the spinous have a spin 
quantum number. Finally, the electron density is de- 
termined via the spinon density since (n-^) = (n®) as 
mentioned earlier. Notice that in the usual slave-boson 
formulation of the tJ model, the doubly occupied state 
is forbidden because U/t — > oo. By comparison, dou- 
ble occupancy is allowed in the slave rotor theory, which 
makes it possible to investigate the effect of charge fluc- 
tuations for arbitrary values of U/t. The slave rotor for- 
mulation was used by Florens and Georges to study the 
Mott transition of the Hubbard model (using single site 



/ J \ -1,(7 



6, 



(8) 
(9) 



where we have defined the doping 6, with S > {S < 0) 
representing hole (electron) doping. For example, the 
electron ground state |^) = l^'e)!^'/), and our task re- 
duces to find the normalized wavefunctions l^t^g), j^*/) 
subject to the constraint that {'$g\n^\'^g) = S and 

We define reduced spinon and rotor Hamiltonian as 



i,j.(T i 

He = (*/|i/sR|*/) 

= -E*^.X..e-^%+^'^^+^E"'("'-l) 

i,j.{T i 

+ ^E(sf-s,^>/- 

Here we have used the notation, (...)/ = (^/| 
and (. . .)g = {"^gl . . . I^e), and also defined 



(10) 



B,, = (e-'^'e+^^O 

Xij = {fiafja)f- 



(11) 
l*/> 

(12) 
(13) 



Note that there is no implicit summation over spin 
a in defining Xij^ and we have assumed Xtj is spin- 
independent. We will denote the nearest and next- 
nearest neighbor values of Bij as B and B' respectively. 
Similarly, the nearest and next-nearest neighbor values 
of Xij will be denoted by x and x'- 

The ground state energy in mean field theory is Eq = 
{^f\Hf\^f) = {^glHgl-^g). Thus, in order to minimize 
Eq, we must choose |^'/) to be the ground state of Hf, 
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and I'i'g) to be the ground state of Hg. This means we 
must self consistently solve for the ground state of the 
two coupled Hamiltonians 

Hf = - Y}^,B,,flf^^ +J J2 Sf • S^-^/ J2 (14) 
He = -2Y}^oX^,e-''^e^''^+^Y.^<f-t^eY]^l (15) 

zj' i ia 

where we have introduced chemical potentials for the 
spinous (/ij) and rotors (//e) as Lagrange multipliers to 
impose the mean field number constraints Eq.® and 
Note that the two decoupled Hamiltonians Hf, Hg are, at 
this stage, still strongly interacting models, and we will 
use further approximations in the spinon and rotor sec- 
tors to make progress. However, note that if we develop 
better techniques to solve these independent spinon and 
rotor Hamiltonians, we can obtain better results within 
this self-consistent mean field approach. 

The ground state energy of the electronic model is 
given by the expectation value of the Hamiltonian i/sR 
of Eq. (O , which leads to 

+ J^(Sf.sJ),. (16) 



IV. SPINON SECTOR 

The spinon sector corresponds to a Hamiltonian with 
kinetic energy and an antiferromagnetic exchange inter- 
action. In principle, the spinon sector can be treated 
using diagrammatic techniques, especially when the spin 
interaction strength is small compared to the spinon ki- 
netic energy. Away from this regime, for larger spin in- 
teractions, the diagrammatic methods may not be ad- 
equate. Here we adopt a simpler mean field approach 
to the spinon sector which has the advantage that we 
can describe the normal state, various broken symmetry 
states, and singlet pairing states of spinous. 



A. Uniform spin liquid and 120° Neel state on the 
triangular lattice 

There have been extensive theoretical stud- 
ies of the Mott transition on the triangular 
lattice2ii2£i22i2^iSai^i^i22i. Among the suggested in- 
sulating states are a uniform spin liquid which can be 
viewed as a 'Fermi liquid' of spinon quasiparticlea22i^, 
and the 120° Neel ordered antiferromagnet. In order to 
study the Mott transition into such states, we decouple 



the exchange interaction as 




+ e^'^-V^^/^^+c-Q-./^^), (17) 

where Xij is real and AI is the sublattice magnetization 
which lies in the XY plane. The wavevector Q is chosen 
to correspond to the 120° Neel state, namely Q • a = 
Q • 6 = 27r/3 where a = x and b = —x/2 + y\fij2 are 
unit basis vectors on the triangular lattice. If M — 0, the 
ansatz reduces to a Fermi sea of spinous, while a nonzero 
M leads to an antiferromagnetically ordered state. 

B. D-wave pairing ansatz 

In order to describe the superconducting state of the 
cuprates, we make a mean field ansatz for the spinon 
Hamiltonian, Hf, following Kotliar and Liu—. We decou- 
ple the antiferromagnetic exchange term in the particle- 
particle and well as in the particlc-hole channel. Upto 
constant terms, this yields 

Here we assume Xij is real, and |Aij| = is mag- 

nitude of the d-wave pairing, with Ajj being positive on 
x-bonds and negative on y-bonds of the square lattice. 
This decoupling scheme highlights the fact that there is 
no difference, at short distance, between the tendency 
to antiferromagnetism and the tendency to form d-wave 
pairs: both involve spin-singlet formation. The difference 
between these states only becomes clear on longer length 
scales, and will be dealt with in a subsequent paper. 

C. Commensurate spin-density wave ansatz 

As a candidate for the underlying normal state of the 
underdoped cuprates, we examine a strongly correlated 
spin-density wave (SDW) state with ordering wavevector 
Q = (tt, tt), which is consistent with that of the adjacent 
antiferromagnetic phase. To describe the commensurate 
SDW ordered at Q, we decouple the exchange interaction 
as 

+ (-1)-+^- y (nj^ - n^^^ - + . (19) 

Here M{—M) refers to the expectation value of — 
(n| — ?^x)/2 on the A{B) sublattices. 
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D. U(l) staggered flux or DDW ansatz 

Another candidate normal state for the weakly doped 
cuprates is the so-called U(l) staggered flusi^i^ or DDW 
state^i^. This state has orbital currents on the bonds 
of the square lattice which form a staggered pattern and 
break time-reversal and lattice symmetries. The mean 
field ansatz for this state is obtained by decoupling the 
exchange interaction as 

Sf •Sj^-^(x, jU- + h.c.). (20) 

where Xji — Xij is complex, with a phase corresponding 
to an enclosed flux of < <&/ < tt per plaquette. <&/ is 
staggered (alternating in sign) from one square plaque- 
tte to the next in a checkerboard pattern. Further, we 
assume Xij is real on diagonal bonds which connect sites 
of the same sublattice. 



V. ROTOR SECTOR: CLUSTER MEAN FIELD 
THEORY 

The rotor Hamiltonian can be solved by a variety of 
techniques used commonly for bosonic systems. Here 
we formulate a self-consistent cluster mean field theory 
which is a straightforward extension of the single-site 
mean field theory used to describe the Mott transition 
in the Bose-Hubbard modeL-" ! - ^ . The idea of the cluster 
mean field theory is to focus on a finite cluster of sites 
and treat the influence of the sites outside the cluster (the 
"bath" ) using a mean field order parameter. The hop- 
ping terms that couple sites within the cluster to sites 
outside the cluster are decoupled using 

e-^e.g*e, ^ e"'*'$, i e cluster, j e bath, (21) 

where the mean field superfluid order parameter $ = 
(e*^') induces number fluctuations in the cluster, and has 
to be determined self-consistently. This is schematically 
shown in Fig. [l]for the square lattice. 

In the superfluid phase, $ ^ 0, and we choose it to 
be real. The Mott insulator phase, on the other hand, is 
characterized by $ = and a charge gap. All terms 
within the cluster, including the onsite U and intra- 
cluster hopping, are retained completely. This cluster 
mean field Hamiltonian is diagonalized exactly to obtain 
the eigen-energies and the ground state rotor wave func- 
tion. The explicit cluster Hamiltonian is given for specific 
cases below. We note that by going to larger clusters, 
we can obtain a progressively better description of the 
bosonic excitations of the rotor sector. 
Single-site mean field theory: Applying the mean 
field theory to a single-site "cluster" , Hg for a square lat- 
tice model is reduced to a single site rotor Hamiltonian, 

= -8itx+t'x'me-^'^+^<^)+^{nr~f^on'. (22) 




FIG. 1: Single site (left) and two-site cluster (right) for the 
suqare lattice. Solid line shows nearest hopping within the 
cluster. Dashed lines represent coupling to the bath described 
by the order parameter $. 



For given x,x' (obtained from the spinon problem) and 
charge density (n^) = 5, we start from some trial value 
of $ and diagonalize H^^^ in the rotor number (angular 
momentum) basis {I^t.^)} by truncating the rotor Hilbert 
space to < L. L is increased until the results con- 
verge. We found L = 2 is accurate enough for all of our 
calculations. Using the ground state wavefunction |G)(i), 
the ground state average and $ = (e*^)(i) are 

computed, fie is adjusted to give the correct charge den- 
sity, and the procedure is iterated until <& is converged. 
Alternatively, one may minimize the ground state energy 
under the charge density constraint {n^){i) to find $, and 
the result will be the same. Once $ is determined self- 
consistently, we record energy eigenvalues Ea and eigen 
states \a). Within the single site approximation, the ro- 
tor kinetic energy B = B' = These values will be 
fed into the spinon Hamiltonian as parameters. For the 
triangular lattice model with t' = 0, the corresponding 
single-site Hamiltonian is 

ijf = -mx^e'^' + e+^') + ^{ny - y^gu' . (23) 

Two-site cluster: Going to larger cluster size yields 
much better approximations for the rotor kinetic energy 
and intersite correlations. For example, in a two site clus- 
ter model of the the square lattice problem, the Hamil- 
tonian Bg is approximated by 

^ _2ix(e-^'^e+'''^+h.c.) 

- (6tx + 8i'x2)$(c"''^^ + c-*^^ + h.c.) 

+ + + (24) 

Similar to the single site case, this Hamiltonian is diag- 
onalized numerically in basis {|7^i,?^2)} ^'^'i *^ is deter- 
mined self-consistently via $ = |(e='=*^i)(2)| = |(e^'^^)(2)|, 
where the average is with respect to the ground state 
|G)(2). Intersite correlations are obtained from the 
ground state wavefunction using B = {e~'^^^e~^'^^^)(2) for 
nearest neighbors and B' = $^ for next nearest neigh- 
bors. These will be used in the spinon Hamiltonian. For 
the triangular lattice, the two site cluster Hamiltonian is 
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given by 



A, (2) 



-2<x(e-'''^e+*''^ +h.c.) 



10tx$(e 



h.c.) 



+ ji<r + j{ntf^MnUn',). (25) 

Finally, we note that the rotor hopping txij is complex 
in the U(l) staggered flux state. The rotor Hamiltonian 
in this case can be solved similarly on a two-site lattice 
to obtain the complex Bij . It is straightforward to gen- 
eralize this to larger clusters depending on the available 
computational resources. 

The main qualitative difference in the results obtained 
from the single-site mean field theory and the cluster 
mean field theories is the following. In single site mean 
field theory, the kinetic energy B is the same as the 
square of the order parameter, Thus, it vanishes 

in the Mott insulating phase. This amounts to neglect- 
ing all number fluctuations within the Mott phase, and 
it is too crude an approximation especially when close 
to the Mott transition. By contrast, the cluster mean 
field theory does much better job capturing the short 
distance correlations, the kinetic energy oc B does not 
simply follow <I>^. For instance at the commensurate fill- 
ing (nf) — 1 and large U/t, deep in the Mott insulator, 
we find B ^t/U. 



VI. MOTT TRANSITION ON THE 
TRIANGULAR LATTICE AT HALF-FILLING 

A Mott transition refers to an interaction induced tran- 
sition from a conducting phase to an insulating phase. 
The insulating phases may or may not be accompanied 
by any conventional broken symmetries. One interest- 
ing recent example of such a Mott transition is the pres- 
sure induced insulator to metal transition in organic com- 
pounds such as K-(ET)2Cu2(CN)3 at an effective density 
of one electron per site (half-filling) The simplest 

model Hamiltonian describing this material is the Hub- 
bard model on an isotropic triangular lattice. Here we 
use the cluster slave rotor mean field theory to study the 
T = phase diagram of model Eq. ([T]) on a triangular 
lattice with t' = 0. 

Previous exact diagonalization studies of the triangu- 
lar lattice Hubbard model suggested^^ a direct first order 
transition from a uniform metallic phase into an antifer- 
romagnetic insulator with 120° Neel order at U/t ~ 12. 
However that study also pointed out^^ that a more com- 
plex phase diagram emerges from Hartree-Fock calcula- 
tions which show a metallic SDW and incommensurate 
spiral ordering in addition to the non-magnetic metal and 
the antiferromagnetic insulator. However, the validity of 
such a weak-coupling approach is unclear for such large 
values of U/t. More recent numerical studies^^ using the 
so-called "path integral renormalization group" (PIRG) 
method indicate that the ground state first undergoes 



a metal to spin-liquid insulator transition at a critical 
Uc ~ 5t. At larger U > Uc, the spin liquid phase un- 
dergoes a second transition into the 120° Neel ordered 
antiferromagnetic state. Variational wave function stud- 
ies, which attack the problem from the insulating phase, 
find that the 120° Neel ordered insulator is unstable 
to a spin liquid state which may resemble a projected 
Fermi sea of spinous^. This spin liquid state with a 
Fermi sea of spinous has been argue d^^i'^^ to provide a 
description of the non-magnetic insulating phase seen in 

K-(ET)2CU2(CN)3. 

Using the ansatz Eq. p?)) for the triangular lattice, we 
can obtain the mean field spinon Hamiltonian 



'k(T 2 
k,o- k 



= —2{tB + ^^)(cosfcQ -i- coskb + COS fee) — ^/ 



(26) 

/• 

(27) 



where ka = k • d, fcf, = k • 6, and fcc = ka + kt- The 
spinon kinetic energy has contributions from the bare 
spinon hopping as well as from the antiferromagnetic in- 
teractions. 

The excitation energy of the SDW quasiparticles is 
given by 



-^k - — 



k-Q 



^^k= ^(ek-Ck-Q)2 + (3JAf)2 



(28) 
(29) 



The self-consistency conditions for M, x and the number 
constraint equation (n/) = 1 — S (which determines /i/) 
reduce to 

^ = ^E^[0M.)-Q(-^^)]' (30) 

k 

^ = ^E["ke(-A,) + .^e(-A+)] 



X (cos ka + COS kf, + cos kc) , 



(31) 
(32) 



Here Q{x) is the step function, and the coherence factor 



1 



Ck - Ck-Q 



1 



(33) 



The results obtained for the Mott transition with a 
fixed J = t/A are plotted in Figl2l We show results for 
the rotor order parameter $, which distinguishes metallic 
($ 7^ 0) from insulating ($ = 0) states, and the sublat- 
tice magnetization M for the 120° Neel state. Using a 
single-site cluster for the rotors, and computing the en- 
ergy of the different states, we find the phase diagram 
in Figl^^a) which shows a direct, strongly discontinuous. 
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FIG. 2: The rotor order parameter <E> and the sublattice mag- 
netization M as function of U /t for the triangular tU J model 
at half filling, J — t/4. Left (a): Predictions of the single-site 
slave rotor mean field theory. The transition from the uni- 
form metal to the antiferromagnetic insulator (AFI) is first 
order. Right (b): Results obtained using a 2-site rotor clus- 
ter. There is an intermediate spin liquid state in which M is 
zero within our numerical accuracy. 



transition from the uniform metal into a magnetically or- 
dered insulating state, occurring at U/t ^ 6. The mag- 
netic ordering in the insulator is 'classical' in the sense 
that the ordered moment is the full moment (M = 1/2), 
and the electronic kinetic energy vanishes in the insula- 
tor. By contrast, the 2-site cluster for the rotor yields 
a richer phase diagram as shown in Figl^fb). It exhibits 
a continuous Mott transition into an insulating state at 
about U/t ^ 7.8. However the insulating state has no or- 
dered moment, M is essentially zero within our numerical 
accuracy. The insulating phase is thus consistent with a 
featureless 'spin liquid' described by a spinon Fermi sea. 
This 'spin liquid' state may be further unstable to spinon 
pairing but we have not, thus far, explored such a pos- 
sibility. For larger U/t ~ 18.5, we find a discontinuous 
transition of this uniform 'spin liquid' state into a 120° 
ordered antifcrromagnet. This result, of a transition at 
large U/t from the spin liquid into the ordered state, can 
be checked using PIRG calculations extended to larger 
values of U/t than in Ref. '^t'. Our main message is that 
using a cluster approach to the rotor sector leads to a 
significant qualitative difference with respect to a single- 
site mean field theory in that it leads to an insulating 
spin liquid phase at intermediate U/t. 

This phase diagram is not too sensitive to the precise 
value of J/t so long as it is not too large. For smaller 
values of J/t, the window of the 'spin liquid' phase ex- 
tends to larger U/t. For much larger values of J/t, we 
find that even the larger clusters can yield a strong dis- 
continuous transition from the metal to a magnetically 
ordered insulator. The details of the phase diagram in 
the slave rotor mean field theory are determined by the 
competition between the renormalized spinon kinetic en- 
ergy, ~ t^ /U , and the fixed exchange constant J. If the 
Mott transition of the rotors occurs at a small U/t, then 
the renormalized spinon kinetic energy is large compared 
to the antiferromagnetic exchange constant J and this 
can frustrate magnetic ordering. This is similar in spirit, 
although not in details, to local charge fluctuations frus- 
trating magnetic order on the triangular lattice via ring- 



exchange processes. Using a variational approach to the 
ring-exchange model, Motrunich has shown that it can 
significantly suppress or completely destroy magnetic or- 
der leading to a uniform spin liquid stat&~. A phase 
diagram for the triangular lattice Mott insulator with an 
intervening spin liquid phase is also consistent with re- 
sults from numerical approaches such as the PIRG. How- 
ever the metal-insulator transition point obtained in that 
study is at U/t ^ 5 indicating that our mean field result 
may not be quantitatively accurate. 



VII. D-WAVE SUPERCONDUCTIVITY IN THE 
DOPED CUPRATE MATERIALS 

High temperature superconductivity in the hole doped 
cuprates has been a prototypical example of how strong 
correlations lead to a variety of novel phenomena not 
observed in conventional superconductors. There are 
two, seemingly quite different, approaches to d-wave su- 
perconductivity in the cuprates. One approach starts 
from an electronic Fermi liquid and shows that anti- 
ferromagnetic spin fluctuations centered at momentum 
(tt, tt) on a square lattice can mediate d-wave pairing of 
electrons^!^i2S,. The other point of view, proposed ini- 
tially by Anderson,— begins from the undoped insulat- 
ing state of the cuprates and argues that this is close to 
forming a spin liquid with d-wave singlet pairing which 
naturally leads to d-wave superconductivity upon doping 
away from half-filling. 

From the perspective of the coupled spinon-rotor the- 
ory these two approaches may be reconciled at low en- 
ergy. The weak coupling spin-fluctuation calculations 
would, when applied to the spinon sector of the slave ro- 
tor mean held theory, predict an instability of spinous to- 
wards d-wave pairing arising from antiferromagnetic fluc- 
tuations induced by J. Thus, although conventional spin 
fluctuation theories and models such as the spin-fermion 
modeP^ are quite different from spin liquid based ap- 
proaches, they appear to be formally identical when rein- 
terpreted as applying to spinous rather than electrons. 
At low dopings, when the exchange interaction becomes 
much stronger than the spinon kinetic energy, the weak 
coupling 'spinon pair' may smoothly cross over into the 
resonating valence bond singlets of the spin liquid Mott 
state. Naively, this crossover is expected to happen when 
the kinetic energy of the doped carriers becomes compa- 
rable to the exchange energy between neighboring spins. 
Recent cellular DMFT calculations^ suggest that this 
crossover may exhibit some features of 'local quantum 
criticality'. 

There has been extensive work on the tJ model using 
U(l) and SU(2) slave boson mean held (and correspond- 
ing gauge) theories^iiSiiii. The slave boson mean field 
theory successfully captures some essential features of the 
high temperature superconductors, including the d-wave 
symmetry, the pseudogap energy scale, and the super- 
conducting dome"'^^. Despite the seemingly bold approxi- 
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FIG. 3: Doping dependence of d-wave superconductivity aris- 
ing from doping a Mott insulator. Left: Maximum supercon- 
ducting gap Aq. Right: Superconducting order parameter 
(CiTCj,i>- 



FIG. 4: Quasiparticle properties of the d-wave superconduc- 
tor. Left: the quasiparticle weight Z obtained from slave- 
rotor mean field theory. Right: nodal Fermi velocity Vf as 
function of doping. 



mations involved in such a slave particle formulation, the 
qualitative predictions of the slave boson mean field the- 
ory agree remarkably well with more sophisticated meth- 
ods like the variational wave function approacb^i^. 

In this section, we apply the cluster slave rotor mean 
field theory to study the superconducting phase of the 
hole-doped square lattice tU J model. Our first purpose 
is to check that it can qualitatively reproduce the well 
known results of the tJ model in the regime U ^ J. 
Second, we show that by solving the rotor problem, we 
learn about the incoherent spectrum weight associated 
with charge fluctuations which leads to, among other 
things, particle-hole tunneling asymmetry. 

The d-wave pairing ansatz for the spinon sector leads 
to the well known Hartree-Fock-Bogoliubov Hamiltonian 

= E^k/L4.-E(^k/^T/-ki+h.c.),(34) 

k,cr k 

= -2(iB + 3Jx/4)(cosfc^ + cosfey) 
— it' B' cos kx cos ky ~ ^ f , (35) 
Ak = Ao(cos /cj; — cos fcy)/2. (36) 

H'^^'^ is diagonalizcd by the standard Bogoliubov trans- 
formation to yield quasiparticle excitation energies = 
V^k + ^k- Here Ao = 3J| Ay I is the antinodal gap. Ak, 
X, and fif are determined self-consistently by 

k " 

1 = -^{cos kx- COS kyf, (38) 

A: 

k 

After self-consistency is achieved, 

X'- (27V)-i^(l-a/Sk)cosfc,cosA:j, (40) 

k 

is computed as input to the rotor sector. 



We have numerically solved the coupled rotor-spinon 
problem for U = 12t, J = t/3, t' = — i/4 (these values 
are typical for hole doped cuprate high temperature su- 
perconductors) at doping < 5 < 0.4 using a two site 
rotor cluster. Fig. [3] shows the doping dependence of the 
antinodal spinon gap Ag as well as the superconduct- 
ing order parameter {ci-\Cj,i) = $^Ao/3J. The spinon 
gap decreases monotonically with increasing doping and 
corresponds to the psuedogap observed in high temper- 
ature superconductors. The doping dependence of the 
order parameter, on the other hand, takes a dome shape, 
reaching maximum at optimal doping 5 ~ 0.2. This is 
in accordance with the doping evolution of Tc of high 
temperature superconductors. The demise of supercon- 
ductivity at low doping is due to the vanishing of con- 
densate density, which in the slave rotor formulation is 
proportional to $^ and goes roughly as 5 at low dopings. 

A. Doping dependence of quasiparticle weight and 
nodal Fermi velocity 

Interaction effects result in a transfer of spectral weight 
from sharply defined coherent quasiparticles to incoher- 
ent excitations. This reduces the quasiparticle weight Z 
from unity. Within the slave rotor mean field theory, Z 
is given by and plotted in the left panel of Fig. ID 
Our slave rotor result is close to the slave boson result, 
as expected since U /t ^ 1 and the number (charge) fluc- 
tuation is suppressed. In the right panel of Fig. 2] we 
plot the nodal quasiparticle Fermi velocity, vp = d^k/dk 
for k ^ kp along the nodal direction. Here in order to 
measure vp in unit of eV-A we have used t = 350meV 
and a lattice spacing 3. 8 A. We find that vp is a, smooth 
function of S and does not vanish as (5 — > 0. The doping 
dependence of Z and vp is in qualitative agreement with 
experiments. 

However, slave boson and slave rotor theories un- 
derestimate Z,vp when compared with Gutzwiller pro- 
jected variational wavefunction calculations^i*^, indi- 
cating that gauge fluctuations that bind slave particles 
have to be taken into account to describe those calcula- 
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tions. Furthermore, the calculated value of vp is smaller 
(by a factor of two) compared with experiments^i^, 
while the doping dependence oi vp in the range S ~ 
0.1 — 0.2 disagrees with the experimental observation^ 
that Vp — 1.5eV-A is essentially completely doping in- 
dependent in this range. On the other hand good agree- 
ment with experiments is achieved by variational wave 
function calculations^i*^. These disagreements show a 
shared weakness of the slave-particle mean field theories. 
Namely some degree of spinon-chargon interaction be- 
yond mean field theory, arising from gauge fluctuations, 
needs to be retained in order to make a quantitative com- 
parison with the photoemission experiments. 



function 



B. Particle- hole asymmetry in tunneling 



Tunneling experiments in the underdoped cuprate su- 
perconductors find a marked asymmetry between the 
tunneling conductance at positive bias (tunneling in an 
electron from tip to sample) versus negative bias (ex- 
tracting an electron from the sample). This had been 
observed in earlier experiments in slightly underdoped 
materials^i^ which focussed on the tunneling conduc- 
tance at low energy. More recent experiments, which 
have gone to lower doping and to larger bias, have shown 
that this asymmetry grows significantly in the strongly 
underdoped regimei^. 

This asymmetry was explained by calculations using 
slave boson mean field theory, which pointed out that it 
arises from the proximity of the underdoped cuprates to 
a correlated insulating state at zero doping^^. Rantner 
and Wen^iS calculated the expected tunneling spectrum 
assuming that the slave bosons behave as free particles 
with a bandwidth ^ t. Here we use the self-consistent 
slave rotor mean field theory to calculate the tunneling 
spectrum into the weakly doped superconducting state of 
the cuprates. We show that the tunneling spectrum has 
a coherent BCS-like form at low energies while it has sig- 
nificant incoherent weight at negative bias in agreement 
with earlier results. We also compute the ratio of the 
integrated (upto a cutoff) tunneling density of states at 
positive to that at negative bias and compare the doping 
and cutoff dependence of this ratio with recent tunneling 
experiments 
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Within the slave rotor mean field approach, using Eq. 
(2]) and (O, the electron Green function factorizes into 
spinon and rotor Green functions. This leads to a sim- 
ple qualitative picture for the electronic tunneling pro- 
cess — namely, tunneling an electron involves tunneling 
a spinon and a charge. This factorization preserves the 
sum rule on the electron spectral function A(k, oj) that 
dujA(k,uj) = 1. The resulting local electron Green 



/3 



Ge(r,zc^„) = - / dTe'^-^c^JT)cU0)) (41) 
Jo 

= - j|''dre-"-l^(4^(r)/t,(0))(e^«-Me-^«-W). 



In another word, on the mean field level, in order to tun- 
nel in an electron one needs to tunnel in both a spinon 
and a charge. Using the Hartree-Fock-Bogoliubov mean 
field solution of the spinous and cluster mean field solu- 
tion of the rotors, we find Ge(v,iajn oj + iS) at zero 
temperature: 



1 

N 



E 



|(m|e-*^^|0)|2 



'^u + i6 - E'k 

|2 



I e 



m 



'u + i6 + Ei,+ (€r, 



eo) 



(42) 



where |m) refers to the m-th eigenstates of the rotor clus- 
ter with eigenenergy em, with |0) being the ground state, 
r is a site belonging to the cluster, = (1 + Ck/-Bk)/2 
and = ^ — are the usual BCS coherence factors. 
This leads to the electronic tunneling density of states 
(TDOS) N{r,Lj) = -7r-i3Gr,c(w), 



N{r,uj) 



1 [ul\{m\e-^'''mH{u-E^-emo) 

5{Lj+E^+e,no)\ ■ (43) 



N 



2|/™l„+j 



where Cmo = Em — co- The first term proportional to 
<i>^, called the coherent part, simply reflects the Bogoli- 
ubov quasiparticle spectrum with energy i^k- It describes 
processes in which the tunneling charge joins the ground 
state (the condensate) with no energy cost. This leads to 
the usual V-shaped nodal spectrum at low energy char- 
acteristic of the d-wave symmetry. The second term is 
the incoherent part, corresponding to processes creating 
an incoherent charge excitation whose local energetics is 
well captured by the cluster mean field theory, and it 
dominates at high energies. 

The calculated TDOS is shown in the left panel of Fig. 
[5]for U = 12t, J = i/3, and doping x = 0.06. Within the 
cluster mean field approximation, the incoherent weight 
consists of discrete peaks corresponding to the excitation 
energy levels of the rotor cluster. We have included a 
small Lorentzian broadening for these discrete levels. We 
find that the tunneling spectrum exhibits a pronounced 
particle-hole asymmetry: there are more significant inco- 
herent charge excitations at negative bias than large pos- 
itive bias. Roughly speaking, upon extracting an electron 
from the system, the nearby electrons which are within 
a few lattice spacings can lower their kinetic energy, by 
~ t, by rattling around the created empty space. On the 
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FIG. 5: (a): The tunneling density of states (TDOS) at dop- 
ing X = 0.06. (b): The dependence of the particle-hole asym- 
metry ratio R (defined in Eq. I44|) on the cutoff energy Qc for 
a series of dopings. 



pendent). This 'partial sum rule' can then be used to 
estimate the hole density locally from tunneling exper- 
iments on doped Mott insulatorsiSiffi. The slave boson 
and rotor mean field theories violate this 'partial sum 
rule' since the spinons and rotors are assumed to be in- 
dependent particles at mean field level. For U/t — > oo, 
the ratio R thus saturates, as a function of fJc, to a value 
which is close to 6{1 + 6) / {1 — S)'^ . This violation is ex- 
pected to be fixed by including gauge fluctuations about 
the mean field theorj*^. 



VIII. CORRELATED NORMAL STATES IN A 
WEAKLY DOPED MOTT INSULATOR ON THE 
SQUARE LATTICE 



contrary, it costs a lot of energy to add electrons to the 
system unless one finds an unoccupied site, and hence 
one has little incoherent spectral weight on the positive 
low bias side. Qualitatively similar ideas have been pro- 
posed by others using variational wavefunctions or exact 
diagonalization of the t-J model on small chistcrs'^^. 

This particle-hole asymmetry can be quantified by in- 
troducing the ratio^ of the positive to the negative 
spectral weights at a given space point r, 



i?(r) 



J^^N{r,E)dE 
f_^N{v,E)dE' 



(44) 



where the cutoff ilc ^ U . 

The right panel of Fig. O shows the fie dependence 
of the ratio R in the doping range x = 0.06 — 0.2. For 
small ric, it starts off close to unity for all dopings, re- 
flecting the nearly particle-hole symmetric d-wave TDOS 
at low bias. Very rapidly, however, the incoherent weight 
present on the occupied side leads to a suppression of R 
which plateaus at a value smaller than unity (~ 0.7) with 
little doping dependence. This is in rough agreement 
with the data in Ref. 18,. At bias 1 < eV/t < 3 there is 
a rapid change in R versus the cutoff, beyond which all 
the curves appear to saturate. This is consistent with the 
bulk of the incoherent weight on the occupied side being 
centered at an energy ~ 2t. (The wiggles in R versus VLc 
in FigIS]are due to the finite size cluster method.) 

It has been shown out that for electronic states which 
obey a Gutzwiller projection constraint (two electrons are 
forbidden from occupying the same lattice site), there is 
a rigorous 'partial sum rule'— which constrains the spec- 
tral weight of occupied and unoccupied states when there 
is a separation of scales t <^ Vic ^ U . In the context of 
tunneling measurements, this translates^ into the result 
that the ratio of the integrated TDOS at positive bias to 
that at negative bias is 2(5/(1 — 5) where 5 is the density 
of doped holes^i^. Furthermore, the ratio i?(r) has the 
advantage that the unknown tunneling matrix elements 
that determine the relation between the TDOS and the 
measured conductance, which can depend on the coor- 
dinate r, cancel out^ (assuming they are energy inde- 



The finite temperature normal state of underdoped 
high temperature superconductors, for T > Tc, is the 
least understood part of the phase diagram. It exhibits 
pseudogap phenomenon and highly anomalous thermo- 
dynamic and transport properties. It has been suggested 
that understanding the underlying T = normal state in 
this regime might shed light on the pseudogap phenom- 
ena. Very recent experiments by Doiron-Leyraud et 
have begun to access this ground state in clean crystals 
of underdoped oxygen-ordered YBa2Cu306.5. Applying 
high magnetic fields, in the range ~ 40-60T, to suppress 
superconductivity, they discovered quantum oscillations 
in the electrical resistance at low temperatures. This dis- 
covery can be most simply understood by postulating a 
well-defined Fermi surface made of small pockets in the 
field induced normal ground state. The authors argue 
that the data is less consistent with a specific band struc- 
ture feature^ of ortho-II YBCO, and more in line with a 
state with broken symmetry giving rise to Fermi surface 
reconstruction. This state is supposed to compete with 
superconductivity, to be close in energy to the supercon- 
ducting state, and to show up as the ground state when 
superconductivity is strongly suppressed. 

Several such candidate states have been proposed. One 
is the commensurate SDW at ordering vector (tTjTt), 
which has been argued to be relevant to the electron 
doped cuprate superconductors^i. Given the proximity 
to the antiferromagnetic Neel phase at low doping, it is 
likely that SDW is stabilized and extends to higher dop- 
ing under high fields. This is supported for example by 
the observation of antiferromagnetism inside the vortex 
cores^i^ and from theoretical studies of the vortex struc- 
ture using Bogoliubov-deGennes (BdG) equationsSI. An- 
other candidate is the U(l) staggered flux statei^, (also 
known as the DDW state^i^). This state is also well 
known to be close in energy with the d-wave supercon- 
ducting stated and has been advocated by Patrick Lee 
and collaborators to describe the pseudogap phaseii. 

In this section we study the SDW state and the U(l) 
staggered flux state using cluster slave rotor mean field 
theory. We then compute the expected cyclotron mass 
TO* for each state as a function of doping, and compare 
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them with the experiment of Doiron-Leyraud et a^. We 
show that by studying the doping dependence of the m* 
one can, in principle, unambiguously differentiate be- 
tween these two scenarios for the normal state. We also 
suggest a possible resolution for the observed difference''^ 
between the hole density as inferred from the quantum 
oscillation measurements and from the Tr value. 



A. Spin density wave state 



Applying the SDW ansatz Eq. p9|) . we obtain the 
mean field spinon Hamiltonian 



kcr k 

- flik+Ql)^ (45) 
where we have defined 



= -2(tB + 3Jx/4)(cosfca; + cosfcy) 

— it' B' cos kx cos ky — fi f , 

-^k = 2 ^ ^k+q) ± 2^k, 

Ei, - x/(ek-a+Q)' + (4Jm)2. 



(46) 
(47) 

(48) 



The mean field self-consistency conditions then read 

k ^ 



1-^ = ^E[Q(-^k) + 0(-A^)] 



(50) 



1 sr^ fCk- Ck+Q 



—T 



{coskx + cos ky) 



[e(-Ak)-e(-A^ 



(51) 



After self-consistency is achieved, we can calculate the 
second neighbor spinon hopping amplitude using 



X 



^ cos k, cos ky [e(-Ak ) + e(-A+)] . (52) 



and use this in solving the rotor model. 

In Figini we plot the sublattice magnetization M as 
function of doping predicted by the cluster slave ro- 
tor mean field theory. M starts at the classical Neel 
state value 1/2 at half filling and drops monotonically 
with increasing doping and vanishes at S ~ 0.17. The 
doped SDW state has four Fermi pockets centered at 
the nodal points k„ = (±7r/2, ±7r/2). Naively, from a 
weak-coupling picture, one expects additional electron- 
like pockets centered around the (tt, 0) point for small 
doping. We find that these do not appear in our self- 
consistent ground state, suggesting that the weak cou- 
pling argument does not work in this regime. Consistent 
with this, the area enclosed by the hole pockets around 
(7r/2, 7r/2) satisfies the Luttinger's theorem for a broken 
symmetry state with a doubled unit cell. 
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FIG. 6: Left; The sublattice magnetization M in the SDW 
normal state. Right: The spinon and rotor fluxes in the U(l) 
staggered flux (or DDW) state. Inset shows the electron cur- 
rent on a nearest neighbor bond in units of et/h versus doping. 



B. Staggered Flux state 

The U(l) staggered fiux ansatz Eq. ([^D]) leads to the 
mean field spinon Hamiltonian 



Ha. 



T^(/k,<T + /k+Q,<T)(/k.'^ ~ h+Q,<r) + h.C. 



(t'B'ji + M/) 5](/L/k,. + /,^+Q,./k+Q..) 



(53) 



where 7^ — 4cosfcj; cos ky, and the complex hopping 



S—:Lx,:Ly 



(54) 



^^fiux is easily diagonalized to yield the the energy spec- 
trum, = -t'B'-fl. -fif±\Dk\. The "+" and the "-" 
branch touches each other at the nodal points k„ re- 
gardless the values of staggered flux, and the low energy 
excitations are massless Dirac fermions. In the doping 
regime of our interest only the "— " branch is occupied. 
The spinon chemical potential /i/ and the complex "or- 
der parameter" x = Xi,i±x = Xli±y are found by solving 
the following self-consistency equations 



i-'^ = i^E[0(-^k) + 0(-^^k)]. 



(55) 



-E 

27V ^ 



D 



cosk,-^[ei-E^) + ei-E+)]. (56) 



Here the sum over k is within the reduced Brillouin zone, 
and N is the number of sites on each sublattice. Finally, 
x' is obtained by 

^'^ ^Eco«^-^o^^yrQ(-^k) + Q(-^k)]- (57) 

k 

With self-consistency achieved in the coupled rotor- 
spinon problem, we obtain x — Ixk*^^'* and B = 
iBle*--/^, where is the spinon flux per plaquette and 
$r is the rotor flux (we work in a symmetric gauge, e.g.. 
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Bi^i±x = B*^^~). The orbital bond current 
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-|xB|sin(^ + -^) 



(58) 

can be taken as the order parameter of the U(l) stag- 
gered flux (or DDW) phase. Figl6] shows the doping de- 
pendence of <&/, <i>r, and the magnitude of the orbital 
current. At half filling the system is in the tt flux state, 
<i>^ and $r have the same magnitude tt but opposite sign, 
accordingly the bond current is zero. The current / in- 
creases with doping, reaching maximum at (5 ~ 0.12. 
For a single plaquette, the maximal bond current cor- 
responds approximately to a magnetic field of 150 Gauss 
at the center of the plaquette. / drops to zero around 
S ~ 0.2, where the U(l) staggered flux state gives way 
to the Fermi liquid state and both and vanish. 
The total area enclosed by the hole pockets increase lin- 
early with doping at small doping and is consistent with 
Luttinger's theorem for a broken symmetry state with a 
doubled unit cell. 



C. Cyclotron mass 

We next discuss whether the quantum oscillation mea- 
surements of Ref. can shed more light on the non- 
superconducting normal state at low doping, and, in par- 
ticular, distinguish between a commensurate SDW metal 
and the U(l) staggered flux (or DDW) state. The quan- 
tum oscillation data offer two pieces of information if the 
field induced normal state is a conventional Fermi liq- 
uid, (although more exotic possibilities cannot be ruled 
ouf*'5.47). The first is the area of the hole pocket which is 
obtained from the oscillation frequency, and the second is 
the cyclotron mass which is obtained from the tempera- 
ture dependence of the oscillation amplitude and which is 
related to the density of states at the Fermi level. In our 
calculations, both the SDW and the U(l) staggered flux 
state have hole pocket areas consistent with the Luttinger 
theorem for a broken symmetry state with a doubled unit 
cell and this information cannot be used to discriminate 
between them. We therefore turn to a calculation of the 
cyclotron mass in these two states to compare with ex- 
periments. 

The cyclotron effective mass m* is related to the den- 
sity of states at the Fermi level in a Fermi liquid state. 



TO* = No 



1^ 



(59) 



Here iVo is the density of states per hole pocket per spin 
at the Fermi level measured in units of \/ta?. As be- 
fore, we set the hopping t = 350meV and the Cu-Cu 
lattice spacing a = 3.8A. The doping dependence of m* 
in the SDW state is shown in the left panel of Fig[71 The 
SDW state is characterized by a diverging to* as 5 ^ 0. 
This reflects the classical character of the magnetically 
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FIG. 7: Cyclotron mass in units of the bare electron mass in 
the SDW state (left) and the U(l) staggered flux (or DDW) 
state (right). 



ordered insulator at 5 = 0, namely the staggered mag- 
netization of the insulator in the mean field theory is 
M — 0.5 and the spinon kinetic energy is zero. This 
leads to a divergent effective mass for 5 — > 0. Spin fluc- 
tuations beyond mean field theory are expected to cut off 
this divergence; nevertheless the mass is expected to be a 
strongly increasing function as 5 decreases. The peak in 
TO* around (5 ~ 0.17 is a feature associated with the van- 
ishing of the SDW order and arises from the large density 
of states around the (tt, 0) point. The magnitude of to* 
in the range 5 = 0.1-0.15 is around 4.5too (toq being the 
bare electron mass), in apparent disagreement with the 
experimental value I.Qtoq. 

The doping evolution of the cyclotron effective mass 
TO* in the U(l) staggered flux (DDW) state is shown in 
the right panel of Fig[7l This state is characterized by 
a square root increase of to*, i.e. to* oc ^/5, at small 
doping. This is a direct consequence of the linear Dirac 
spectrum at low energies and therefore a linear density 
of states. Again, the peak feature at (5 ~ 0.18 is associ- 
ated with the sudden disappearance of the staggered flux 
state. The value of to* ~ Stoq at 5 = 0.1 is closer to the 
experiment result I.Otoq. 

While a comparison of the cyclotron effective mass 
with experiments may seem to favor the staggered flux 
state over the SDW state as a normal state candidate, we 
cannot rule out the SDW scenario based on this quanti- 
tative comparison. As we have seen in previous sections, 
slave particle theory is not quantitatively accurate at the 
mean field level and underestimates the nodal vp in the 
superconducting state by a factor of two at (5 = 0.1. To 
distinguish these two scenarios, it is more useful to study 
the doping dependence of to* . From Fig|7]we observe that 
the effective cyclotron mass in the SDW state and the 
staggered flux state have almost opposite doping depen- 
dences, we propose that measurements of to* at a series 
of nearby dopings, especially at smaller 5, may help to 
identify which candidate normal state is realized. 

Finally, the experiments-y" find a difference between the 
carrier densities as inferred from the quantum oscillation 
period and from T^. Assuming 2 hole pockets in a re- 
duced Brillouin zone, which is consistent with a SDW or 
a U(l) staggered fiux state, the oscillation period leads 
to an estimate 5 ~ 0.075 while the value estimated from 
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Tc is closer to 6 ^ 0.1. This discrepancy^, an appar- 
ently significant violation of the Luttinger count, can be 
accounted for in several ways including a possible bilayer 
splitting^-. However, the bilayer splitting is expected to 
be very small at low doping due to renormalization by 
strong correlation effects. We offer another plausible ex- 
planation for this discrepancy. The field induced 'normal 
state' is likely to contain vortices, which can be inferred 
from the sign change, with temperature, of the Hall co- 
efficient, a feature which is seen in the mixed state in 
other cuprate materials and at other dopings^. It has 
been suggested that the underlying normal state of the 
cuprates should appear in the vicinity of such vortices^. 
However, it has also been found in BdG calculations^ of 
the structure of a single vortex that the density around 
the vortex is not uniform. At low doping and even in the 
presence of moderate Coulomb interactions^, the den- 
sity appears to be depleted within the vortex core and 
piled up as a screening cloud at some distance away from 
the core, with the far region being superconducting and 
at the bulk density. If the normal state being explored 
arises from the depleted core region, the observed den- 
sity could deviate from the bulk density in the correct 
direction. Testing this explanation calls for further ex- 
periments as well as detailed theoretical calculations of 
the vortex structure. 



D. Implications for Fermi arcs 

We end this section by commenting on the relation 
between the hole pockets revealed in the low tempera- 
ture quantum oscillation experiments and the Fermi arcs 
observed in the high temperature normal state above 
Tc in photoemission experiments. Assuming that the 
field induced normal state has hole pockets, it has been 
suggested that the Fermi arcs seen in photoemission 
studies^ may be just these hole pockets. However, the 
quantum oscillation experiments are carried out at rather 
high magnetic fields B ~ 50T at very low temperature 
T <^ Tc, while the photoemission experiments are done 
in zero magnetic field and at fairly high temperatures 
T ^ 40 — 200K . Given this, there is no a priori rea- 
son to expect that the photoemission experiments are 
accessing the properties of such an underlying normal 
ground state. The Fermi arcs seen in the photoemission 
experiments seem to be more consistent with fluctuat- 
ing d-wave superconductivity^!^. In particular, earlier 
work by us^S, suggests an explanation which relies on phe- 
nomenological spin charge separation ideas, quite close to 
the route explored in this paper, and provides a possible 
explanation for Fermi arcs whose length scales T/T*, 
where T* is the pseudogap temperature. Current pho- 
toemission experiments do not appear to have the reso- 
lution to settle the issue of whether the arcs indeed arise 
from hole pockets. We suggest that further high resolu- 
tion photoemission experiments could focus on the mo- 
mentum distribution of the low energy spectral weight 



(integrated from zero to some small energy). It would be 
useful to check whether this spectral weight tends to be 
centered at momenta corresponding to the low tempera- 
ture d-wave nodes of the SC, or whether it tends to shift 
away towards some underlying hole pocket Fermi surface 
once the system is heated to temperatures above Tc- 



IX. DISCUSSION 

We have formulated a cluster slave rotor mean field 
theory of strongly interacting electronic systems, which 
extends earlier work of Florens and Georges -"^^^ In this 
formulation, the rotor sector and spinon sector are cou- 
pled to each other by self-consistency equations, and 
the rotor Hamiltonian is solved on a finite cluster self- 
consistently coupled to an order parameter bath (rest of 
the lattice). The main advantage of the cluster mean 
field theory lies in that the short range correlation func- 
tions of the rotors is properly taken into accomit. By 
applying it to several examples, we have shown that this 
method provides a unified framework to study strongly 
correlated electronic systems. Although we only focused 
on zero temperature, it is straight forward to generalize 
the theory to finite temperature. Extensions to inhomo- 
geneous states are also likely to be of interest. 

One shortcoming of the cluster mean field theory for 
the rotor sector is that, by working with finite clus- 
ters coupled to a static bath, we freeze the Goldstone 
modes (phase fiuctuations). Such Goldstone mode fluctu- 
ations can be explicitly retained in complementary meth- 
ods such as the random phase approximation^''' or other 
approaches^ to the rotor sector. These methods of solv- 
ing the rotor problem can be combined with the spinon 
sector to formulate the self-consistent slave rotor mean 
field theory™. 

In the path integral formulation of the slave rotor 
approach^^'^^'^®, our mean field theory is equivalent to a 
saddle point approximation which ignores the gauge field 
fiuctuations. Such fluctuations beyond mean field theory 
are expected to introduce effective interactions between 
slave particles to yield the correct low energy sum rules^ 
including the 'partial sum rule' for the particle-hole 
asymmetry ratio R. These will be explored elsewhere^. 

Note added: After our paper was submitted, some 
related works have appeared on the arxiv. The self- 
consistent cluster mean field theory for bosons has been 
shown to successfully describe the supersolid phase and 
the phase diagram of bosons on a triangular lattice^^. 
Two further experiments carried out on YBa2Cu408 
(hole doping ^ 0.15) have also found evidence of 
Shubnikov-de Haas oscillations2L2^ in a magnetic field 
with a larger cyclotron mass, m* ~ 2.7 — SrTig, than in 
ortho-II YBCO. A recent preprint by Chen et al discusses 
the field induced antiferromagnetic metallic ground state 
in the underdoped cuprates^*. Their principal result is 
the relation between the quantum oscillation frequency 
and the hole doping in this state and it is in agreement 
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with our conclusion. Various proposals for interpreting 
the Fermi arcs above T^, as well as their connection to 
the hole pockets in ground state, have been examined 
recently by Harrison et al''^ and Norman et al2&. 
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